Derivation of effective parameters of magnetic metamaterials composed of passive 

resonant LC inclusions 
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We determine rigorously the effective permeability of magnetic metamaterials composed of passive 
particles exhibiting magnetic resonances. The effective permeability is expressed in terms of circuit 
parameters characteristic to these particles (L, C, and R) and particle geometry. The derivation 
takes into account the magnetic coupling between particles. We show that a bigger concentration 
of particles leads to an improvement in the metamaterial performances (such as bandwidth or loss 
tangent), but this effect saturates as the particles are packed more tightly. The theory is validated 
through numerical simulations of physically realizable metamaterials. 
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Since the development of the first electromagnetic 
metamaterials based on split-ring-resonators (SRRs) and 
wires, it has been acknowledged that the narrow band- 
width, high dispersion, and significant loss inside them 
have been the main factors that limited their practical 
use in commercial applications. Theoretical analysis of 
the individual particles that cornpose such metamaterials 
have been performed before 0, and gave valuable in- 
sight on their physics. However, a complete and reliable 
theoretical derivation of the effective material parameters 
of these metamaterials in terms of their geometry has not 
been developed yet. 

In this paper we focus on metamaterials having nega- 
tive permeability, which are commonly designed and built 
using periodic arrays of unit cells composed of inclusions 
such as SRRs , omega particles (^] , and other types of 
similar resonant particles 0. Our purpose is to derive 
rigorously the effective permeability of such metamateri- 
als in terms of the circuit parameters characteristic to the 
composing particles (i.e. inductance, L, capacitance, C, 
and resistance, R) and the particle geometry represented 
by the geometry factor F. We will show that larger F 
factors lead to an increase in the performance in terms 
of bandwidth or loss tangent. However, we will show 
that F cannot be made arbitrarily large and reaches an 
upper bound as the particles are packed closer together. 
Numerical simulations will be performed to validate the 
theoretical results. 

Consider the unit cell depicted in Figure ^ AH the 
resonant particles commonly used to generate negative 
permeability (such as SRRs, omega particles) are essen- 
tially LC resonant circuits [l|, Q that can be represented 
in general by a planar, capacitively loaded loop as de- 
picted in Figure n Losses are taken into account by us- 
ing a resistor, R, in series with the capacitor, C, and the 
loop. For now, the loop can have an arbitrary shape. A 
metamaterial is generated by using periodic arrays of this 
unit cell. Consider a plane wave propagating through this 
metamaterial in the +z direction, such that the magnetic 
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FIG. 1: Unit cell containing a capacitively loaded loop. Peri- 
odic arrangements of this unit cell generate the metamaterial. 
Inset: The capacitively loaded loop used to generate the mag- 
netic response of the metamaterial. The curve Cioop follows 
the contour of the loop; Aioop is the area delimited by Cioop 



field, Bi = yBoe^^^"'^ (e^'^* time dependence assumed), 
is perpendicular on the plane of the loop. Writing Fara- 
day's Law in integral form for an arbitrary loop inside 
the metamaterial we obtain: 

f E • dl = -jiu I (Bi + Bioop + 51 B") • (1) 

where Bioop is the magnetic field produced by the cur- 
rent induced in the loop by the incident Bi, Bn is the 
field produced by the current induced in the n-th loop. 
From the definition of auto inductance, J^^ Bioop 'ds = 
LioopI, where L/oop is the auto-inductance of the loop, 
and / is the current induced in the loop. Similarly, 
f. Bn • ds — Mnin, whcrc Mn and /„ are the mu- 

^loop 

tual inductance between the n-th loop and the loop for 
which Eq. is written, and, respectively, the current 
through the n-th loop. Note that M„ can be positive 
or negative depending on the relative position between 
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the two loops. In a bulk material, far from the edges, 
because of the symmetry of the metamaterial structure, 
we can consider that / = /„ for every n. Moreover, the 
voltage across the resistor in series with the capacitor 
is I[R+ 1/ijLuC)] = /p^^^ Bi • ds. Therefore, Eq. ((TJ 
becomes: 
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Bi • ds + IL. 



eq 



(2) 



where Lg// = i + X^n ^^n- Letting ujq = (CLeff)^^^^ be 
the resonant frequency of the loop, the above equation 
gives the current induced in the loop: 



Leff 



jujR/Leff 



(3) 



In general, the magnetic flux through the loop, 
/^^ Bi • ds, is a function of the geometry of the loop. 
However, it is well known that in order for the meta- 
material to be described in terms of effective material 
properties, the unit cell has to be much smaller than the 
wavelength. Therefore, for frequencies around loq where 
the permeability of the metamaterial becomes negative, 
we can make the approximation J^^ Bi-ds ~ BiAioop = 
^oHiAioop- Note that Bi and Hi represent the local in- 
cident B and H fields related through Bi — noHi since 
the loops are placed in vacuum. For example, for square 
loops (i.e. a geometry widely used for building SRRs) 
/^^^^^ Bi • ds = BtAioopSinc{kolioop/'2), where koop is the 
length of the loop. If we require X/koop = 10, as it is 
generally desired, then 1 > sinc(fco/ioop/2) > 0.99, which 
allows us to neglect the sine term. Under this approxima- 
tion, it follows from (PJ that the magnetic dipole moment 
per unit volume is given by 
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' LeffVuc - cjg - jujR/Leff 



(4) 



where Vuc is the volume of the unit cell. 

Note that Mjup- follows a purely lorentzian model, 
with the quality factor of the material given by 



Q = R/iL,ffLUo) = 1/{looRC) 



(5) 



For weak magnetically coupled loops, L^ff ~ L, and Q is 
approximately the quality factor of the individual loop. 
Also, note that the term: 
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LeffVuc 



(6) 



only depends on the geometry of the unit cell and that of 
the loop, and, as recognized in Ref. has an important 
contribution to the performances of the metamaterial. 

The relative effective permeability of the metamaterial, 
given by = 1 + M/Hi^ can be expressed in terms of 
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FIG. 2: The effective permeability retrieved from the S pa- 
rameters (solid line) , and the permeability computed using 
(dashed line). This excellent match was obtained for F = 0.16 
and Q = 71. 



the geometry factor, F, quality factor, Q, and resonant 
frequency, ujq, as: 



Fu^ 



(7) 



To verify the accuracy of (jSJ through Q we performed 
numerical simulations in Ansoft HFSS similar to those 
described in Refs. |2 and0. The metamaterial has one 
cubic unit cell in the direction of propagation, and has 
infinite extent in the transverse direction (this is simu- 
lated using perfect electric and perfect magnetic bound- 
ary conditions). The unit cell has dimensions 40 by 40 
by 40 mm. The square loop has dimensions 30 by 30 
mm and is loaded with a capacitor C = 1.108 pF and 
a resistor R = 3.617. These dimensions give a resonant 
frequency slightly below 500 MHz, which means that the 
unit cell diameter is less then A/15, and the metamate- 
rial can be described by effective material parameters. 
The effective material parameters are retrieved from the 
S parameters using a standard method widely used for 
these metamaterials j^], and are compared with the pa- 
rameters computed theoretically from Eq. Written 
in terms of capacitance, which is assumed to be constant 
(the capacitive coupling between adjacent loops can be 
neglected), and resonant frequency (which can be eas- 
ily retrieved from the S parameters), the geometry fac- 
tor is (from Eq. ©) F ^aCulAj^^^/Vuc- We obtain 
Wo = 27r • 487.5 MHz, and consequently the expected F 
and Q factors are F = 0.17, and Q = 77. The real Q and 
F values can be retrieved by matching the permeability 
determined from the S parameters with that given by 
Eq. (fTI) . As can be seen in Figure [21 an excellent match 
was obtained for F = 0.16 and Q ~ 71. Note that these 
values are within 10% of the expected F and Q, which 
validates the theoretical relations lO through (71). 

Equation |(7J) gives the lorentzian variation expected 
for these magnetic metamaterials Ij. The relative per- 
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meability can be controlled to some extent by controlling 
the three parameters that enters Q. Thus, is rela- 
tively easy to control by loading the loop with lumped 
capacitors as suggested in Ref. 6. The quality factor, Q, 
directly responsible for the losses inside the metamate- 
rial, is harder to control since it depends on the dielectric 
and ohmic losses Out of the three parameters, F is 
the easiest to control since it depends only on the geom- 
etry of the unit cell. Moreover, a closer look to Eq. 
reveals that a bigger F factor leads to an increase in the 
negative permeability bandwidth and a decrease in the 
magnetic loss tangent, in agreement with Ref. Q. Conse- 
quently, the next step is to determine how large F can be 
made. Note that the parameters that enter Eq. © that 
gives F are dependent on each other. The area of the 
loop is strongly related to its inductance, and the unit 
cell volume. If we reduce the unit cell volume by pack- 
ing the loops closer together in an attempt to increase 
F, the effective inductance of the loops would increase 
as a result of an increase between the mutual inductance 
between the loops. This suggests that F cannot be made 
arbitrarily large. 

Since it is commonly used in literature, we will focus 
next on metamaterials made of cubic unit cells. For the 
commonly used square loop A^oop/Kc = koopikoop/lucf , 
where Uoop is the diameter of the loop, and luc is the di- 
ameter of the unit cell. Obviously, the ratio hoop/luc is 
always less than unity, and typically is less than 0.8 to 
reduce the coupling between adjacent loops. For exam- 
ple, in Ref. 10 this ratio is 0.79, while in Ref. ^ it is as 
low as 0.52. On the other hand the inductance is pro- 
portional to ^lokoop- Since the unit cell is cubic we can 
neglect the inductive coupling between loops. Under this 
approximation, a more exact equation written for square 
loops made of cylindrical wires gives 

Lioop = f 2.303 logio 5^ - 2.853) (8) 

Zn \ w J 

where w is the diameter of the wire. Typically, the term 
hoop I w is of the order of hundreds, and has a moderate 
influence on the inductance Ljoop because the logarithm 
is a slowly varying function for large arguments. As- 
suming hoop I w > 100 it follows that Lioop > 3.33^ohoop, 
which combined to the constraint on Af^^^/Vuc gives (see 
Eq. ®) 

F « = 0.15 (9) 

3.33/io/ loop 

The above value was obtained for cubic unit cells and 
square loops. We expect loops of different shapes to give 
slightly lower values of F because the ratio Af^^^/Vuc is 
smaller than for square loops while the inductance re- 
mains proportional to /io\A^oop- For example, for cir- 
cular loops it can be shown in a similar manner that 
F fa 0.12. It follows from the above analysis that for 
cubic unit cells and arbitrary loop shapes F < 0.2. 



This bound on F is determined directly from Eq. ((Hjl, 
and is supported by both the parameters retrieved from 
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FIG. 3: The geometry factor vs the concentration of particles 



the HFSS simulation (shown in Figure|21), and the results 
obtained by others using numerical methods j^. Note 
that this analytically derived upper bound on F is sig- 
nificantly smaller than that predicted by other analytical 
models [l|. 

These results were obtained for cubic unit cells which 
are commonly used in the design of magnetic metama- 
terials. To increase F over the limit determined above, 
we can decrease the unit cell volume by reducing the cell 
axial dimension, or, equivalently, by increasing the con- 
centration of loops. However, besides a decrease in Vuc, 
we also expect a substantial increase in the effective in- 
ductance, Leff, due to the magnetic coupling between 
the loops (recall that F oc l/Lg//). The F factor can, 
thus, be accurately computed either analytically, from © 
by computing the mutual inductances that give ie//, or 
numerically using Ansoft HFSS and the procedure de- 
scribed above. For simplicity, we used the latter method 
to find F versus the concentration of loops (expressed 
as the number of loops per cubic unit cell) . The results 
are presented in Figure 13 and show the geometry factor 
asymptotically approaching an upper bound. Thus, even 
for tightly packed loops (i.e. 20 loops per cubic unit cell), 
F remains below 0.7. 

In conclusion, we derived the analytical equations that 
allow an accurate prediction of the effective permeability 
of magnetic metamaterials composed of passive resonant 
inclusions commonly used in literature, such as the split- 
ring-resonator, or the omega particle. The analysis took 
into consideration the magnetic coupling between these 
particles. We showed that, for the commonly used cu- 
bic unit cell, the geometry factor F is typically 0.15. An 
increase in the F factor can be achieved by increasing 
the concentration of particles. However, even for tightly 
packed particles F reaches an upper bound significantly 
lower than 1. We validated these results through numer- 
ical simulations performed in Ansoft HFSS. 
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